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Abstract 

In light of the recent LHC Higgs search data, we investigate the pair production of a SM-like 
Higgs boson around 125 GeV in the MSSM and NMSSM. We first scan the parameter space of each 
model by considering various experimental constraints, and then calculate the Higgs pair production 
rate in the allowed parameter space. We find that in most cases the dominant contribution to 
the Higgs pair production comes from the gluon fusion process and the production rate can be 
greatly enhanced, maximally 10 times larger than the SM prediction (even for a TeV-scale stop 
the production rate can still be enhanced by a factor of 1.3). We also calculate the x? value with 
the current Higgs data and find that in the most favored parameter region the production rate 
is enhanced by a factor of 1.45 in the MSSM, while in the NMSSM the production rate can be 


enhanced or suppressed (osusy /osy varies from 0.7 to 2.4). 


PACS numbers: 14.80.Da,14.80.Ly,12.60.Jv 


I. INTRODUCTION 


Based on the combined data collected at the center-of-mass energies of 7 TeV and 8 TeV, 
the experimental programme to probe the mechanism of electroweak symmetry breaking 
at the LHC has recently witnessed the discovery of a new particle around 125 GeV n 
The properties of this particle, according to the updated analyses of the ATLAS and CMS 
collaborations at the end of 2012 bi. roughly agree with the the Standard Model (SM) 
prediction and thus it should play a role in both the symmetry breaking and the mass 
generation. However, the issue of whether this particle is the SM Higgs boson is still open, 
and indeed there are some motivations, such as the gauge hierarchy problem and the excess in 
the di-photon channel over the SM prediction Al H to consider new physics interpretation 
of this boson. Studies in this direction have been performed intensively in low energy 
supersymmetry (SUSY) and it was found that some SUSY models can naturally provide a 
125 GeV Higgs boson Lii and fit the data better than the SM |6] (similar studies have also 
been performed in some non-SUS Y models like the little Higgs models and two-Higgs-doublet 
or Higgs-triplet models [7]). 

After the discovery of the Higgs boson, the next important task for the LHC is to test the 
property of this Higgs boson by measuring all the possible production and decay channels 
with high luminosity. Among the production channels, the Higgs pair production is a rare 
process at the LHC. Since it can play an important role for testing the Higgs self-couplings 
li (the determination of the Higgs self-couplings is of great importance since it is in- 
dispensable to reconstruct the Higgs potential), it will be measured at the LHC with high 
luminosity. 

In the SM the Higgs pair production at the LHC proceeds by the parton process gg > hh 
through the heavy quark induced box diagrams and also through the production of an off- 
shell Higgs which subsequently splits into two on-shell Higgs bosons min 
rate is rather low for vs = 14TeV, about 20 fb at leading order [10] and reaching roughly 35 


The production 


fb after including the next-to-leading order QCD correction |8|. The capability of the LHC 
to detect this production process was investigated in 15]. These analyses showed that for 
a 125 GeV Higgs boson the most efficient channel is gg > hh — bbyy with 6 signal events 
over 14 background events expected for 600 fb integrated luminosity after considering 


some elaborate cuts (the detection through other channels like hh > bbW*+W~ and 


hh > bbrt7- has also been studied recently hd hd). In principle, the capability can be 
further improved if the recently developed jet substructure technique [16] is applied for the 
Higgs tagging. 

The Higgs pair production at the LHC may also be a sensitive probe for new physics. In 
supersymmetric models such as the Minimal Supersymmetric Standard Model (MSSM) [17], 
the pair production of the SM-like Higgs boson receives additional contributions from the 
loops of the third generation squarks and also from the parton process bb > H; > hh with 
H; denoting a CP-even non-standard Higgs boson ha mt It was found that in some cases 
(e.g., a light stop with a large trilinear soft breaking parameter A; and/or a large tan 8 
together with moderately light H;), these new contributions may be far dominant over the 
SM contribution, and as a result, the rate of the pair production may be enhanced by several 
orders "m Note that since the experimental constraints (direct or indirect) on the SUSY 
parameter space have been becoming more and more stringent, the previous MSSM results 
should be updated by considering the latest constraints. This is one aim of this work. To 


be specific, we will consider the following new constraints: 


e The currently measured Higgs boson mass m, = 125 GeV B. In SUSY this mass is 
sensitive to radiative correction and thus the third generation squark sector has been 


tightly limited. 


e The LHC search for the third generation squarks " 5o far although the relevant 
bounds are rather weak and usually hypothesis-dependent, it becomes more and more 


clear that a stop lighter than about 200 GeV is strongly disfavored. 


e The observation of B, > utu” by the LHCb . In the MSSM it is well known that 
the branching ratio of B, > ut; is proportional to tan? 8 /m/, for a large tan 8 and 
a moderately light H " Since the experimental value of B, — utpu coincides well 


with the SM prediction, tan 8 as a function of my has been upper bounded. 


e The LHC search for a non-standard Higgs boson H through the process pp > H > 
T'T- . Such the search relies on the enhanced Hbb coupling and the nought signal 
seen by the LHC experiments implies that a broad region in the tan 8 — mg plane has 


been ruled out. 


e The global fit of the SUSY predictions on various Higgs signals to the Higgs data re- 
ported by the ATLAS and CMS collaborations b4. the dark matter relic density 
as well as the XENON2012 dark matter search results can also limit SUSY pa- 


rameters in a complex way. 


Another motivation of this work comes from the fact that the Next-to- Minimal Supersym- 
metric Standard Model (NMSSM) is found to be more favored by the Higgs data and 
the fine-tuning argument |6]. So far the studies on the Higgs pair production in the NMSSM 
are still absent. So it is necessary to extend the study to the NMSSM. 

This paper is organized as follows. In Sec. II we briefly introduce the features of the Higgs 
sector in the MSSM and NMSSM. Then in Sec. III we present our results for the Higgs pair 
production in both models. Some intuitive understandings on the results are also presented. 


Finally, we summarize our conclusions in Sec. IV. 


II. HIGGS SECTOR IN MSSM AND NMSSM 


As the most economical realization of SUSY in particle physics, the MSSM has 
been intensively studied. However, since this model suffers from some problems such as the 
unnaturalness of u parameter, it is well motivated to go beyond this minimal framework. 
Among the extensions of the MSSM, the NMSSM as the simplest extension by singlet 
field has been paid much attention. The differences between the two models come from 


their superpotentials and soft-breaking terms, which are given by 


Wussm = Y,Q- HU — Y4Q- H4D — Y.L - HaÊ + uÉ,- Ha, (1) 
Wanssu = VQ HU — VQ HaD — VL Hal + MÅ Had + 508%, (2) 
VA = mH? + maHa[ + (Bul, : Hat hc), (3) 
VINSSM — må lH. + saa Hal? + ISP + (A ASH, - Ha+ Beng? +h.c.). (4) 


Here Å; (i wd) and $ denote gauge doublet and singlet Higgs superfields respectively, 
Ô, U : D, L and É represent matter superfields with Y; (i — u,d,e) being their Yukawa 
coupling coefficients, Mm; (i = u,d, S), B, Ax, and A, are all soft-breaking parameters and 
the dimensionless parameters Å and « reflect coupling strengthes of Higgs self interactions. 


Note the p-term in the MSSM is replaced by Higgs self interactions in the NMSSM, so when 


the singlet field $ develops a vacuum expectation value s, an effective u is generated by 
Heff — As. 
Like the general treatment of the multiple-Higgs theory, one can write the Higgs fields in 


the NMSSM as 
H+ Va + Qatipa 1 
Ha i ; ; Ha = i v2 D S = s + — (o + i£), (5) 
Vy + Puten Hy V2 


and diagonalize their mass matrices to get Higgs mass eigenstates: 


Hi u E u 
d Ar H7 
H, = Ug Qa , Å» = UA Qd D = Uc i (6) 
Gr Ay 
Hz» Oo G E 


Here Hi, Hə, H3 with convention my, < my, < mg, and A, A» with convention ma, < 
Ma, denote the physical CP-even and CP-odd Higgs bosons respectively, G? and G* are 
Goldstone bosons eaten by Z and W bosons respectively, and H* is the physical charged 
Higgs boson. The Higgs sector in the MSSM can be treated in a similar way except that it 
predicts only two physical CP-even states and one physical CP-odd state, and consequently, 
the rotation matrices Uy and Ua are reduced to 2 x 2 matrices. 

One distinct feature of the MSSM is that Hı usually acts as the SM-like Higgs boson 
(denoted by h hereafter) and its mass is upper bounded by mz at tree level. Obviously, to 
coincide with the LHC discovery of a 125 GeV boson, large radiative correction to ma, is 
needed, which in turn usually requires the trilinear soft breaking parameter A, to be large. 


For example, in the case of large ma and moderate tan B, Mp is given by 


må = M3 cos? 28 + 2 mad (1- = ) : (7) 
where the first term is the tree-level mass and the last two terms are the dominant corrections 
from the top-stop sector, m; = VM m, (mi, denotes stop mass with convention m;, < m;,) 
represents the average stop mass scale and X, = A, — u cot B. One can easily check that for 
a 500 GeV and 1 TeV stop, |A;| should be respectively about 1.8 TeV and 3.5 TeV to give 
mn œ~ 125 GeV. 

In the NMSSM, m, exhibits at least two new features "n One is that it gets additional 


contribution at tree level so that m2 „ee = (m? — A?v?) cos? 28 + A^v?, and for A ~ 0.7 and 
h,tree Z 


tan 6 ~ 1, Mp can reach 125 GeV even without the radiative correction. The other feature is 
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FIG. 1: Feynman diagrams for the pair production of the SM-like Higgs boson via gluon fusion in 
the MSSM and NMSSM with Hr denoting a CP-even Higgs (I = 1,2 for the MSSM and I = 1,2,3 
for the NMSSM) and Gi; (i, j = 1,2) for a squark. The diagrams with initial gluons or final Higgs 
bosons interchanged are not shown here. For the quarks and squarks we only consider the third 


generation due to their large Yukawa couplings. 


that the mixing between the doublet and singlet Higgs fields can significantly alter the mass. 
To be more explicit, if the state Hı is h, the mixing is to pull down the mass, while if Hə 
acts as h, the mixing will push up the mass. Another remarkable character of the NMSSM 
is that in the limit of very small \ and « (but keep p fixed), the singlet field decouples from 
the theory so that the phenomenology of the NMSSM reduces to the MSSM. So in order to 
get a Higgs sector significantly different from the MSSM, one should consider a large A. 
Throughout this work, we require 0.50 < Å < 0.7 in our discussion of the NMSSM and 


we consider two scenarios: 
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FIG. 2: Feynman diagrams for the parton process bb — hh in the MSSM and NMSSM. 


e NMSSMI scenario: Hj, acts as the SM-like Higgs boson. For this scenario, the addi- 
tional tree-level contribution to mp is canceled by the mixing effect, and if the mixing 
effect is dominant, the parameters in the stop sector will be tightly limited in order to 


give mn ~ 125 GeV. 


e NMSSM2 scenario: Hə acts as the SM-like Higgs boson. In this scenario, both the 
additional tree-level contribution and the mixing effect can push up the mass. So for 
appropriate values of Å and tan, m; can easily reach 125 GeV even without the 


radiative correction. 


III. CALCULATIONS AND NUMERICAL RESULTS 


In SUSY the pair production of the SM-like Higgs boson proceeds through the gluon 
fusion shown in Fig[I]and the bb annihilation shown in Fig[2] These diagrams indicate that 
the genuine SUSY contribution to the amplitude is of the same perturbation order as the 
SM contribution. So the SUSY prediction on the production rate may significantly deviate 
from the SM result. To ensure the correctness of our calculation, we checked that we can 
reproduce the SM results presented in and the MSSM results in lid. Since the analytic 
expressions are quite lengthy, we do not present here their explicit forms. 

In our numerical calculation we take m; = 173 GeV, m; = 4.2 GeV, mz = 91.0 GeV, 
my = 80.0 GeV and a = 1/128 bal. and use CT10 to generate the parton distribution 
functions with the renormalization scale up and the factorization scale upr chosen to be 2mp. 
The collision energy of the LHC is fixed to be 14 TeV. Then we find that for m; = 125 GeV, 
the production rate in the SM is 18.7 fb for gg — hh and 0.02 fb for bb — hh (the rates 
change very little when mp varies from 123 GeV to 127 GeV). 
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FIG. 3: The scatter plots of the surviving samples, showing osusy/osm versus the SM-like Higgs 
boson mass. The plus ”+ (blue) denote the results with only the gluon fusion contribution, while 


the circles ’o’ (pink) are for the total results. 


For each SUSY model we use the package NMSSM Tools-3.2.0 to scan over the pa- 
rameter space and then select the samples which give a SM-like Higgs boson in the range of 
125 + 2 GeV and also satisfy various experimental constraints, including those listed in Sec- 
tion I. The strategy of our scan is same as in |6] except for three updates. First, since the rare 
decay B, > utu has been recently observed with Br(B, > utu) = 3.2413 x 107? 21), 
we use a double-sided limit 0.8 x 107° < Br(B, 3 utu”) € 6.2 x 107°. Second, for the 
LHC search of the non-standard Higgs boson, we use the latest experimental data bd The 
third one is that we require stops heavier than 200 GeV bd. After the scan, we calculate 
the Higgs pair production rate in the allowed parameter space. We will demonstrate the 
ratio &:susy /osm for each surviving sample. Of course, such a ratio is less sensitive to higher 
order QCD corrections. 

In Fig. [3] we show the normalized production rate as a function of the Higgs boson mass 
for the surviving samples in the MSSM and NMSSM (for the NMSSM we show the results 
for the NMSSM1 and NMSSM2 scenarios defined in Sec.II). This figure shows two common 
features for the three scenarios. One is that the production rate can deviate significantly 
from the SM prediction: in most cases the deviation exceeds 3096 and in some specail cases 
the production rate can be enhanced by one order. The other feature is that for most 
cases the dominant contribution to the pair production comes from the gluon fusion, which 
is reflected by the approximate overlap of ’o’ (pink) with ’+’ (blue). Fig. [3] also exhibits 


some difference between different scenarios. For example, in the MSSM the bb annihilation 


contribution can be dominant for some surviving samples, which, however, never occurs in 
the NMSSM. Another difference is that the NMSSMI tends to predict a larger production 
rate than other scenarios. 

Now we explain some features of the results in Fig. B] First, we investigate the cases of 
the MSSM where the bb annihilation plays the dominant role in the production. We find 
that they are characterized by a moderately large tan 8 (tan 8 ~ 10 so that the Hbb coupling 
is enhanced), a moderately light H (300 GeV < my £ 400 GeV) and a relatively large Hhh 
coupling. While for the NMSSM scenarios, since we are considering large A case, only a 
relatively small tan is allowed so that the H;bb coupling is never enhanced sufficiently B. 
We also scrutinize the characters of the gluon fusion contribution in the MSSM. As the first 
step, we compare the sbottom loop contribution with the stop loop. We find that for the 
surviving samples the former is usually much smaller than the latter. Next we divide the 
amplitude of Fig. [I] into five parts with Mj, M», M3, M4 and Ms denoting the contributions 
from diagrams (1)+(2), (3)+(4), (5), (6)+(7) and (8)+(9)+(10), respectively. For each of 
the amplitude, it is UV finite so we can learn its relative size directly. We find that the 
magnitudes of Mə and M3 are much larger than the others. This can be understood as 
follows: among the diagrams in Fig. [I] only (3), (4) and (5) involve the chiral flipping of the 
internal stop, so in the limit m;,,m;, > 2m; the main parts of M» and M3 can be written 


as 


A? A? 
M ^ o5 Y? (c sin? 20,— + o—) (8) 

Tn z 

ti t2 
where Y, is the top quark Yukawa coupling, 0, and A, are respectively the chiral mixing 
angle and the trilinear soft breaking parameter in the stop sector, and c; and c are O(1) co- 


efficients with opposite signs. Since a large Å; is strongly favored to predict m; ~ 125 GeV 


in the MSSM [3] and the other contributions are usually proportional to m?/m? or m; /mz ; 
one can easily conclude that Mə and Ms; should be most important among the five ampli- 
tudes. In fact, we checked that without the strong cancelation between Mz and Ms, the 
production rate can easily exceed 100 fb for most surviving samples. 

As a proof for the validity of Eq.(8), in Fig. 4] we show A;/m;, versus mz, where the 
samples are classified according to the value of R = esusv(gg > hh)/osu(gg > hh). The 
left panel indicates that in the MSSM the region characterized by a light m;, and a large 


|A,/ m; | usually predicts a large R. This can be understood as follows. In the MSSM with 
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FIG. 4: Same as Fig [3] but showing A;/m;, versus m;,. The samples are classified according to 
the value of R = osusy(gg > hh)/osm(gg — hh) with o denoting the hadronic cross section via 


gg — hh. 


a light t4, the other stop (t2) must be sufficiently heavy in order to predict m; ~ 125 GeV 
3]. Then, after expressing sin? 20, in terms of A; and stop masses, one can find that the first 
term in Eq.(8) scales like (A;/ mi, )*(m?mz / ms ); and therefore its value grows rapidly with 
the increase of |A;/m;,| and is unlikely to be canceled out by the second term in Eq.(8). In 
fact, the upper left region of the panel reflects such a behavior. This panel also indicates 
that even for tı and tz at TeV scale, the production rate in the MSSM may still deviate from 
its SM prediction by more than 30%. This is obvious since |A;| in Eq.(8) is usually larger 
than stop masses i Finally, we note that for m; > 1 TeV, there exist some cases where 
the deviation is small even for A;/m;, ~ 3. We checked that these cases actually correspond 
to a small mass splitting between tı and f». In such a situation, the first term in Eq.(B) is 
proportional to A?/ må (since 0, ~ 7/4), and its contribution to the rate is severely canceled 
by the second term. 

Eq.(8) may also be used to explain the results of the NMSSMI scenario. In this scenario 
we checked that the mixing effect on m, often exceeds the additional tree level contribution 
(as discussed in Sec. II), and consequently the soft breaking parameters in the stop sector 
are more tightly limited than the other two scenarios. For example, given the same values 
of m;, and mz, for the three scenarios, the NMSSMI scenario usually prefers a larger |A|. 
Consequently, this scenario tends to predict the largest production rate according to Eq.(&8]. 
As for the R value in the NMSSM2 scenario, the situation is quite complex because a large A 


alone can push the value of m; up to about 125 GeV and thus the soft breaking parameters 
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FIG. 5: Same as Fig. B| but showing osusy/0sm versus X?. Here only the samples satisfying 125 


GeV x mn < 126 GeV are plotted. 


in the stop sector are not so constrained by the Higgs mass B. But, anyway, this scenario 
still has the features that R is maximized for a large A; and a light £; and that R can deviate 
sizably from unity for TeV-scale stops. 

Finally, we focus on the samples which predict a SM-like Higgs boson in the best fitted 
mass region, 125 GeV < mp < 126 GeV [24] 
with the LHC Higgs data (for details, see 
rate osusy/osm in Fig. [5] This figure indicates that in the MSSM and NMSSM2 scenarios, 


. For these samples, we calculate the x? value 


||) and show its correlation with the normalized 


there exist a lot of samples with x? much smaller than its SM value (x2,, = 16.5), which 
implies that the MSSM and NMSSM2 scenarios may be favored by the current data [6]. In 
contrast, the NMSSMI scenario can only slightly improve the fit. From this figure we also 
see that in the favored parameter space with a small y? the production rate can sizably 
deviate from the SM prediction ( in the parameter space with a large x? the production rate 
can be several times larger than the SM value). For example, in the low x? region of the 
MSSM, the normalized rate is approximately 1.45, while in the NMSSM2 scenario the rate 


varies from 0.7 to 2.4. 


IV. SUMMARY AND CONCLUSIONS 


Recently, the CMS and ATLAS collaborations announced the discovery of a new reso- 
nance whose property is in rough agreement with the SM Higgs boson. But the nature of 


this new state, especially its role in electroweak symmetry breaking, needs to be scrutinized. 
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So the most urgent task for the LHC is to test the property of this Higgs-like boson by 
measuring all the possible production and decay channels with high luminosity. Among the 
production channels, the Higgs pair production is a rare process at the LHC. Since it can 
play an important role for testing the Higgs self-couplings, it will be measured at the LHC 
with high luminosity. 

In this work we studied the pair production of the SM-like Higgs boson in the popular 
SUSY models: the MSSM and NMSSM. To make our study realistic, we first scanned the 
parameter space of each model by considering various experimental constraints. Then we 
examined the Higgs pair production in the allowed parameter space. We found that for 
most cases in both models, the dominant contribution to the pair production comes from 
the gluon fusion process with its rate maximized at a moderately light f, and a large trilinear 
soft breaking parameter A+. The production rate can be sizably enhanced relative to the 
SM prediction: osysy/ogsm can reach 10, and even for a TeV-scale stop it can also exceed 
1.3. For each model we also calculated its X? with current Higgs data and found that in the 
most favored parameter region the value of osysy/osm is approximately 1.45 in the MSSM, 


while in the NMSSM it varies from 0.7 to 2.4. 
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